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A RIGOROUS GEOMETRIC DERIVATION OF THE CHIRAL ANOMALY 

IN CURVED BACKGROUNDS 

CHRISTIAN BAR AND ALEXANDER STROHMAIER 


Abstract. We discuss the chiral anomaly for a Weyl field in a curved background and show 
that a novel index theorem for the Lorentzian Dirac operator can be applied to describe the 
gravitational chiral anomaly. A formula for the total charge generated by the gravitational 
and gauge field background is derived directly in Lorentzian signature and in a mathemat¬ 
ically rigorous manner. It contains a term identical to the integrand in the Atiyah-Singer 
index theorem and another term involving the ^-invariant of the Cauchy hypersurfaces. 


1. Introduction 

Anomalies in quantum field theory appear as a violation of conservation laws of currents, 
in other words, as currents that are classically preserved but whose quantum counterparts 
are not. These anomalies are of direct physical significance. A prominent example is the 
chiral anomaly (also ABJ-anomaly or axial anomaly). It explains the rate of decay of the 
neutral pion into two photons ([1, 8]), 7To —> 77. We would like to refer to the monograph [10] 
for further details and applications in quantum field theory. Apart from high energy physics, 
it has also been proposed that this anomaly can be observed more directly in crystals (see 
[34], and also [39] for recent experimental evidence). 

Some historical remarks. Whereas anomalies were first discovered in perturbative com¬ 
putations in quantum field theory, their appearance is related to the Atiyah-Singer index 
theorem. Indeed, the perturbative computation yields a term that looks precisely like the 
local Chern character form of the connection induced by an external electromagnetic field 
and perturbative computations on curved space-times yield a term that is formally identi¬ 
cal to the Atiyah-Singer integrand ([27, 17]). However, these computations were performed 
in Lorentzian signature whereas the Atiyah-Singer index theorem for Dirac operators is a 
theorem in Euclidean signature. 

A direct computation of the divergence of the anomalous chiral current and the relation to 
the local index theorem was studied in Euclidean signature by Nielsen, Romer, and Schroer 
in [35], and by Dowker in [19] in great detail. Another popular method using the Euclidean 
formulation of quantum field theory and path integrals is the Fujakawa method ([20]), which 
has also been applied in curved backgrounds ([21], see also [2]). Both explanations of the 
relation to index theory can be seen only as formal manipulations because the process of 
Wick rotation, passing from Lorentzian to Euclidean signature, cannot be made sense of in 
general. 
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Another, more mathematically rigorous, way to understand anomalies in quantum field 
theory in external fields is via a careful analysis of the second quantization procedure ([11]). 
This was done by Klaus and Scharf for fermions in an external field ([28, 29]). The one- 
particle time-evolution in an external held can mix the negative and positive energy solutions 
of the Dirac equation. This mixing can occur in such a way that the vacuum is mapped to 
a charged sector under the second quantized time-evolution ([28, 29, 31, 13]). The charge 
generated by this process relates to the Fredholm index of a certain operator constructed out 
of the positive energy projection and the time-evolution operator. This is the point of view 
we will take in this paper. 

The fact that fermions and chiral charge can be created by a gravitational held was 
discovered by Gibbons in [22] (see also [24]). He computed this number in several models and 
showed that a Wick rotated Euclidean version was described by the Atiyah-Patodi-Singer 
index formula [23], i.e. by the Atiyah-Singer integral and a correction term involving the 
^-invariant of the boundary. As pointed out in [22], passing from the Euclidean formulation 
to the Lorentzian is highly problematic as only few Lorentzian manifolds admit an analytic 
continuation to a real Riemannian space which is non-singular. This is also mentioned by 
Lohiyah in [30] where the loss of chiral charge from the Reissner-Nordstrom black hole was 
computed directly but found to be consistent with the predictions of the Atiyah-Singer 
formula. 

The present article. We will show in this paper that the gravitational chiral charge creation 
is given by the index of the Lorentzian Dirac operator with certain boundary conditions. We 
use a novel index theorem for hyperbolic operators proved by the authors in [7] to show 
that this index is indeed, in analogy to the Atiyah-Patodi-Singer index theorem in Euclidean 
signature, given as a sum of the Atiyah-Singer integrand and a correction term involving the 
^-invariant. This therefore establishes the formula directly for Lorentzian space-times and 
avoids the use of a Wick rotation. This way we overcome the above mentioned problems with 
analytic continuation and provide a direct explanation why the Atiyah-Patodi-Singer formula 
yields the correct result in the exactly solvable models. Note that hyperbolic operators such 
as the Lorentzian Dirac operators have significantly different analytic properties from elliptic 
Dirac operators in Euclidean signature. For example, the fact that the hyperbolic Dirac 
operator with spectral boundary conditions is a Fredholm operator and has a well-defined 
index at all, is due to the global propagation of singularities theorem and not to local regularity 
theory as in the elliptic case. 

In order to make a precise mathematical statement about the created physical charge 
we work in the framework of algebraic quantum field theory in curved space-times. This is 
natural as the notion of particle is not covariant and it is more appropriate in this context 
to speak about observables and states. We start by explaining in detail how the Weyl held 
is quantized in a globally hyperbolic space-time with spin structure. As is usual in quantum 
held theory in curved space-times we split this procedure into two parts (see for example [38] 
for an introduction). The held algebra can be constructed in a functorial (covariant) manner 
from the held equation. This was hrst done by Dimock ([18]) for the Dirac held. 

In the second step we look at physical representations of the held algebra induced by 
states. The Hadamard condition singles out a class of states whose behavior resembles that 
of hnite energy states in Minkowski space-time. In our context the important fact is that 
Hadamard states differ only by smooth integral kernels. Quantities that are a priory singular 
and need regularization, such as the expectation value of the energy momentum tensor or 
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the electric current, make sense as relative quantities between two Hadamard states. It is 
therefore possible to define the relative energy-momentum tensor and the relative current 
between two Hadamard states without the need for regularization. 

A particular example of a Hadamard state can be defined if the space-time is ultrastatic in 
a neighborhood of a smooth spacelike Cauchy hypersurface. Near this Cauchy hypersurface 
a one parameter group of time-translations can be defined and it makes sense to define the 
vacuum state with respect to this time-evolution by employing the usual frequency splitting 
procedure. We will think of this state as the vacuum state seen by an observer moving along 
the distinguished timelike Killing vector field near the Cauchy hypersurface. If there are two 
such Cauchy hypersurfaces we can compare the two vacua. Due to the interaction with the 
gravitational background and the external gauge field between the two Cauchy hypersurfaces 
it can happen that the difference of the expectation values of the electric current between 
the two states is non-zero. The difference of the expectation values of the total charge 
operator turns out to be given essentially by the index of the Lorentzian Dirac operator with 
Atiyah-Patodi-Singer boundary conditions. We then apply the index theorem proved by the 
authors in [7] to obtain an explicit formula for the charge generated by the background fields. 

More recent related work. An index formula for the scattering operator of the Dirac 
equation coupled to an external field was established by Matsui in [32, 33]. Bunke and 
Hirschmann [12] generalized this to also allow time-dependent metric connections. In these 
results there is no ^-correction term because they are related to the special case of our result 
where the future and past spacelike Cauchy hypersurfaces coincide so that the ^-contributions 
cancel. 

The computation of the divergence of the anomalous regularized current can also be 
carried out directly in Lorentzian signature using the Hadamard regularization. This method 
makes use of the unique singularity structure of Hadamard states and subtracts the singular 
part from the singular current to make it a well-defined covariant and local object. The 
computation of the divergence of the resulting current was done in [40] for the mixed and 
pure gravitational anomalies and in [15] for the trace anomaly that also can be described by 
an index theorem ([14]). The algebraic expressions coincide with the ones obtained in the 
Euclidean framework and those obtained in perturbative computations. The divergence of 
the anomalous current captures aspects of renormalization in curved space-times. However, 
the actual generation of charge in external fields is described by the index of the Dirac opera¬ 
tor and the r;-correct ion term is vital to guarantee that the total charge generated is an integer. 
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2. The setup 

We start by describing the classical setup which gives rise to the one-particle evolution. 
We will consider massless fermions with internal symmetries in an external held on a curved 
space-time which satisfy a Dirac equation. Let (X,g) be an even-dimensional Lorentzian 
manifold where g has signature (+, We assume that (X, g) is globally hyperbolic 

meaning that there exist Cauchy hypersurfaces so that there will be a well-posed initial value 
problem for the Dirac equation. Then (X, g) is time-orientable and we fix one time-orientation. 
Moreover, we assume that (X,g) is spatially compact, i.e. that the Cauchy hypersurfaces are 
compact. We assume a spin structure on X is given so that the complex spinor bundle 
SX -A X is defined. 

Finally, to model the internal symmetries, a Hermitian vector bundle E — > X is fixed and, 
to describe an external gauge held, a compatible connection V E on E is given. Compatibility 
means that the Leibniz rule holds for X E and the scalar product on E. For example, E 
could be a Hermitian line bundle; then V E is a U(l)-connection and can describe the electro¬ 
magnetic potential. 

2.1. The Dirac operator. Spinors are sections of the bundle SX 0 E. Let V be the con¬ 
nection on SX 0 E induced by the Levi-Civita connection X s on SX and the connection X E 
on E, i.e. Vx(ff 0 e) = (V^cr) 0 e + a 0 Vfe. The Dirac operator acts on spinors and is 
locally given by 

if = ig ap 1 a Vp 

where we used Einstein’s summation convention and the coefficient matrices satisfy {7a > 7 / 3 } = 
2 g a p. Here and henceforth {•, •} denotes the anticommutator. 

Since the dimension of X is even the spinor bundle splits into the subbundles of left-handed 
and right-handed spinors, SX = SlX 0 SrX. The Dirac operator interchanges chirality, i.e. 
with respect to the splitting SX 0£ = SlX 0E0 SrX 0 E we have 



Note that there is no mass term on the diagonal. 

The bundle SX comes equipped with a natural nondegenerate but indefinite inner product 
(•,•). We use the convention that (•,•) is antilinear in the first argument and linear in the 
second. The subbundles S l /rX are isotropic with respect to this inner product. The same 
is of course true for the induced inner products on SX 0 E and S l /rX 0 E. Recall that the 
inner product on E is positive definite. 

The Dirac operator is formally selfadjoint with respect to (•, •), i.e. 

(1) [ (i^u,v)dV = [ (u,ifv)dV 

Jx Jx 

where u, v € Cq°(X; SX 0 E) and dV is the volume element on X induced by the Lorentzian 
metric. 

The formally dual operator of f R : C°°(X; SrX 0 E) -A C°°(X; SlX 0 E) is the operator 
f£ : C°°(X; S* L X 0£*) A C°°(X; S^X 0£*) on the dual bundle characterized by 

(2) f (%v)(u)dV = [ v(f R u )dV 

Jx Jx 

where u € Cfi°(X; SrX 0 E) and v € Cfi°(X; S* L X 0 E*). 
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Example 2.1. If E is a Hermitian line bundle, then the curvature of the dual connection on 
E* is the negative of that of E. Up to a factor i, the curvature is a real 2-form which describes 
the electro-magnetic field. Locally, if we write V® = d a + iA a then Vf* = d a — iA a . For the 
Dirac operators this means i*j7 = i^ a {d a + iA a ) = i^ a d a — A^ 01 while i= Ef a d a + A a 7 °. 
Thus i^7 should be thought of as the operator with charge opposite to that of i'tf. 

2.2. Dirac conjugation. For u € SrX x <g> E x define u G S* L X x (g> E* by u(v) = (u,v). 
The map u i—>• u, SrX <g> E —>• S* L X <S) E* is antilinear. Similarly we get an antilinear map 
7 : SlX ig> E —> SftX <S> E* and the inverse maps S* R j L X <g> E* Sl/rX <g> E are also denoted 
by 7 . Equations (1) and (2) imply 

(3) if*u = (■ if)u . 

Dirac conjugation is defined as the map 

(4) r :SIX®E*®S r X®E^S* l X®E*®SrX®E, u®v^v®u. 

Clearly T is antilinear and an involution, T 2 = 1. 

2.3. The Cauchy problem. Denote the space of smooth solutions of ^rU = 0 by SoI(Vr) := 
{u € C°°(X;SrX <g) E) | rf R u = 0} and similarly for % j . V/), and %*. The Cauchy problem 
for the Dirac equation on globally hyperbolic manifolds is well posed (see e.g. [18, Thm. 2.3]). 
This means that for any smooth spacelike Cauchy hypersurface Eel the restriction map 

px:Sol(%)^C 00 ('E ] S R X®E), u^u\r, 

is an isomorphism of topological vector spaces. If E and Yf are two smooth spacelike Cauchy 
hypersurfaces of X, then we put 

Uiv s := pv o (ps)’ 1 : C°°(S; SrX <g> E) -A C' 00 (E / ; SrX ® E). 

This evolution map Uri s extends to a unitary isomorphism 

(5) Us',5 : L 2 (S; SrX <g> E) -A L 2 (S'; SrX ® E ). 

Hence there is a unique Hilbert space completion SoI(Vr) of SoI(Vr) such that for each 
smooth spacelike Cauchy hypersurface Eel the restriction map extends to a Hilbert space 
isometry 

Pr:So\(%)^L 2 (^SrX®E), 

see [6, Lemma 3.17]. We denote the scalar product on SoI(Vr) by (•,•). Similarly, we 
obtain the Hilbert space completions SoI(Vl), SoI(Vj-j), and SoI(Vl) of SoI(Vl), SoI(Vk), 
and Sol(^*), respectively. 

2.4. The fermionic propagator. Let Gr : Cfi°(X; SlX ® E) — A C°°(X; SrX ® E) be the 
difference between retarded and advanced fundamental solutions of i^7r- Then Gr maps onto 
the space of solutions Sol(^) of the Dirac equation ^’ R u = 0. This operator is sometimes 
called the fermionic propagator of i^7 R . Similarly, we obtain linear maps 

C L : CZ°(X- SrX ® E) - Sol^) C C°°(X- S L X <g> E), 

Gr ,* : C 0 °°(X; S* R X ® E*) -» Sol® C C°°(X; S£X ® E*), 

G l ,* : C 0 °°(X; ® E*) - Sol® C C®; ® E*). 

Equation (3) implies 

(6) 


Gr^u = Gru and Gl,*v = Glv 
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for all u G Cfi°(X-, SlX®E ) and v G C^{X\ SrX®E). Using an integration by parts, we can 
express the scalar product on the Hilbert space SoI^r) as follows (compare [18, Prop. 2.2]) 


(7) ( G R v,f) = - [ (v,f) dV 

J x 

where v G C^°(X]SlX © E) and / G SoI(Vr). Analogous formulas hold for the scalar 
products on Sol(^), Sol(^), and Solffll), respectively. 

By [5, Thm. 4.3] the fermionic propagators extend to operators on distributional sections, 


e.g. 


G r : <V°(X; S L X © E) -A C-°°(X; S R X © E). 

The image is precisely the space of distributional solutions of the Dirac equation. 


2.5. Fermionic propagator and Cauchy problem. Let E C X be a spacelike smooth 
Cauchy hypersurface with future-directed unit normal field ns. Any u G C' _00 (E; SlX © E) 
gives rise to a distribution uds on X via (u6s)(v) = f s (psv)(u) dA for all v G C°°(X ; S* L X © 
E*). Since the support of u5 s is contained in E and hence compact, we can apply Gr to 
u5s- An argument using integration by parts shows that / = Gr(u5 s) G C~°°(X; SrX <g> E) 
is the solution of the Cauchy problem 

(8) i%f = 0 and psf = is u - 

Here i/ij, = n^a denotes Clifford multiplication by ns in accordance with Feynman’s slash 
convention. The restriction Gr : L 2 (T,;SlX © E) C Cq°°{X\ SlX © E) — > Sol(^) is 
therefore the inverse of the Hilbert space isometry i/i-^ps : SoI(Vr) — > L 2 (E;SrA" © E) and 
hence is itself a Hilbert space isometry. 


3. States and the chiral anomaly 

In order to describe the quantized Weyl field we start with the construction of the held 
algebra. It is essentially the CAR algebra associated with the space of solutions of the Dirac 
equation. More precisely, let /C be the Hilbert space sum Sol(^) © So1($r). We denote the 
scalar product by (•, • )/c • This Hilbert space comes with the Dirac conjugation T : /C —>■ 1C as 
defined in (4). Then the selfdual CAR algebra associated with the pair (/C,T) is the unital 
^-algebra generated by symbols B(f ) where / G 1C and relations 

(9) / i-a B(f) is complex linear, 

(10) {B(f),B(g)} = (Tf,g) K , 

(11) B(f)* = B(Tf). 

The CAR algebra admits a unique C*-norm and we will denote its C'*-completion by A. This 
is our held algebra. 

3.1. The field operators. We define the field operators by 

T : C 0 °°(A; S* R X © E*) -A A, (u) = B(Gr,*u 0 0), 

T : C 0 °°(A; S L X © E) -A A, W(v) = B(0 0 Grv). 

Relations (9)—(11) together with (6) and (7) imply that 

u i-a 'l'(ti) is complex linear, 
u G f (r) is complex linear, 
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{^Oi), ^(u 2 )} = o, {'I'Oi), ^(t> 2 )} = o, 

=-i f (v, Gr^u) dV, 

Jx 

<F( u )* = 0(u), 

( 12 ) = 0 , = 0 . 

Here = 0 is to be understood in the distributional sense, i.e. 'Sf(')t' R u) = 0 for all u £ 
Cg°(X; S* L X 0 £*) and similarly ^(%v) = 0 for all v £ C^(X\ S R X 0 £). 

3.2. The n-point functions. A state w on this algebra is of course uniquely determined by 
its n-point functions 

:=w(5(/i) ••■£(/„)). 

This also defines a distributional section w n £ C' _00 (X n ; Kl n (S’^A 0 E 0 S* L X 0 if*)) on the 
n-fold Cartesian product of X by 

0n('IH 0 * * * 0 Un) •— k^n(G?/i, • • • , G tin) 

where we have used the notation G = G Rt * 0 G R . We will also refer to this distribution as 
the n-point function. 

A state is called quasi-free if for all n = 1, 2, 3,... and £ K, 

^2n-l(/l, • • • , /2n-l) = 0, 

n 

n(n — 1 ) ^ ^ ^ ^ 

^2n(/l, • • • ,/2n) = (“1 }^ Sgn ^ 

J=1 

where the sum is taken over all permutations s of {1,... , 2n} such that 

s(l) < s(2) < • • • < s(n), s(j ) < s(j + n). 

This means that the state is completely determined by its two-point function. 

3.3. Hadamard forms and relative currents. A two-point function is said to be of 
Hadamard form if its wavefront set WF(0 2 ) satisfies 

WF(0 2 ) C {(xi,£i,X 2 ,£ 2 ) £ T*(X x X) | (aq,£i) ~ (x 2) —£ 2 ), £2 is future directed}, 

where (xi,£i) ~ (x' 2 , £ 2 ) means that these vectors are in the same orbit of the geodesic flow. 
It is known (see e.g. [37, 26] in the case of the Dirac field) that Hadamard forms are unique 
up to smooth kernels, i.e. if uj\ and uj 2 are states with two-point functions of Hadamard form 
then 

w 1)2 := wi — w 2 £ C°°(X x X] M 2 (S r X 0£® S* L X 0 E*)). 

Note that 0i >2 is a smooth bi-solution of the Dirac equation. We can then define the relative 
current J ^ 1,0J2 £ fl 1 (A) as follows. Given a state uj we define S w : C^°(X] S R X 0 E*) —> 
C~°°(X-, S* L X 0 E*) as the unique operator such that 

u('S’(v)'i>(u)) = (Suu)(v) 

for all v € Cq°(X; SlX 0 E). Note that if uj has Hadamard form then we actually have S u : 
Cq > (X\ S r X®E*) —> C°°(X-, S r X®E*). We put S := S U1 — S U2 . If both oj\ and uj 2 have two- 
point functions of Hadamard form then the operator S has smooth integral kernel. We denote 
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the integral kernel of S evaluated at ( x, y) G MxM by s(x, y) G Horn( S* R X y 0 E*. S* L X x ®E*). 
Moreover, by (12), we obtain 

(13) %o S = S o% = 0. 

For £ G TX x denote Clifford multiplication with £ by £ G Hom(5/jAf x 0 E x . SlX x (g) E x ) and 
its dual by £* G Honr(S'2X a; 0 E*, S* R X x 0 E*). Now £* o s(x, x ) G End(S , ^X x 0 E*) and we 
can set 

J W1 ' U3 (0 ■= tr(jf os(x,x)). 

Since s is smooth on M x M, this defines a smooth one-form J^ 1 ^ 2 g fl 1 (X). In physics 
terminology, one could write this definition as 

= 1™ ^i(^ A (a;)(7 /i )|^ s (y)) -uj 2 {^ A (x)^ ll )f^ B (y)) S j , 

where the Einstein summation convention was used on the dotted spinor index A and spinor 
index B. The relative current can be thought of as the expectation value W 2 O J^(x):) of the 
normally ordered current operator where the normal ordering has been done with 

respect to state lo\. It follows directly from the definitions that 

WJ]=W and [/, f*] = w* 

where / G C°°(X ) is a function and V/ its gradient vector held. For any smooth compactly 
supported function / G Cq°(X) we get, using (13) twice, 

[ fSJ" 1 '" 2 dV = [ ) dV = Tr (W*S) = Tr {f%S - %fS) = -Tr (fS%) = 0. 

Jx Jx 

Here we have used the canonical trace Tr on the algebra of integral operators whose integral 
kernel have compact support in the first variable. Its properties are described in Appendix A. 

This shows 5J U1,U2 =0. If £ is a spacelike smooth Cauchy hypersurface and ns the 
future-directed unit normal vector held along £, then we can integrate the relative current 
J Wl ’^ 2 (ns) along the Cauchy hyper surface. Since the current is co-closed the integral does 
not depend on the choice of Cauchy hypersurface and dehnes the relative right-handed charge 
Q^’ U2 between two states 

(14) Q™ := J J^ 2 (n E (x))dA(x). 

3.4. Fock representations. If "H is a Hilbert space the fermionic Fock space J-('H) is defined 
to be the Hilbert space direct sum 

OO 

E{W> = 0A k n. 

k=0 

As usual the vector 1 G C = A°"H C X(J-i) is denoted by Q and is called the vacuum vector. 
On this space we have the usual creation and annihilation operators 

a\v)h = v A h, 
a(y)h = i v h. 

Note that a^(v) is linear in v, whereas a(v) is antilinear in v. 

Let P : X —> X be an orthogonal projection such that 


p + rpr = 1. 
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This induces an orthogonal splitting 1C = P © YP where P = PIC. We obtain a faithful 
♦-representation n of A on the Fock space P{P) by setting 

:= a) {Pf) + a(PYf). 

This representation is called the Fock representation. One checks that for such a representa¬ 
tion the state 

uj(-) = {Cl, ■ Cl) 

is a quasifree state and the Fock representation is canonically isomorphic to the GNS represen¬ 
tation of uj. See [3] or [36] for a detailed discussion of CAR algebras and Fock representations. 

3.5. Fock representations constructed near S. If T is a spacelike smooth Cauchy hy¬ 
persurface then the restriction map 

ps : 1C = Sol{%) 0 SoI(Vr) ->• /Cs := L 2 (T-S* L X 0 E*) 0 L 2 (S; S R X 0 E) 

is an isomorphism of Hilbert spaces. We can therefore identify the Hilbert space 1C with 
conjugation Y with /Cs and conjugation T|s. Moreover, we can identify >S7 ;X|e with ST,, the 
spinor bundle of T. Using the isomorphism ifi^ : SlX\z —> •S'rA'Ie we also get an identification 
of SlX |s with ST. The Dirac operator Vs on T anticommutes with ifi s and is selfadjoint 
with respect to the positive definite scalar product •) on ST. This implies 

(15) %v = 

In particular, conjugation 7 maps Vs-eigenspaces for positive eigenvalues to Vs-eigenspaces 
for negative eigenvalues and vice versa. We consider the spectral projectors 

P>{%) ■= X[o,oo) (Vs) and p>{%) := X(o,oo)(^e) ■ 

Similarly, we define p<(Vs)j P< (Vs); an d the corresponding projectors for Now 

P:=P>{%)(B P> {%) 

is an orthogonal projection on 1C s and (15) implies P + TPr = 1. Its Fock representation 
will then be modeled on the one-particle Hilbert space 

P^ = P+®P~ 

where 

n+=p>{%)L 2 (T-,S*T®E*), 

^s = P>(%)L 2 {T-, ST®E)= P< (%)Li(T-,S*T®E*). 

The space P^ is the one-particle space for particles and P ^ is the one-particle space for 
antiparticles. The Fock space then splits into a tensor product 

-F('Hs) = E{P^)®F{P^). 

As usual we denote by c and cl the creation and annihilation operators on P{P^) and by d 
and rft the creation and annihilation operators on P(P^). Each of them acts on the entire 
Fock space by extending it by the identity map on the other tensor factor. The representation 
of A in this Fock representation is explicitly given by 

vr s (T(u)) = c* (p>(fx)PzG R ,*u) + d{p > {%)p^G R u), 

7te(^(u)) = c(p>(^)pzG R) *v) + d t (p<(V 7 s)psGiju). 
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The notation 

7TE (£(/)) = a\PpY.f) + a(PTpxf) 
with the substitution / = G(u © v) is however more compact. 

It is well known (see for example [16, Sec. 4]) that the vacuum expectation value wj with 
respect to this representation is a Hadamard state, i.e. its two-point function is of Hadamard 
form, if the metric on X and V E are of product type near £. This means that a neighborhood 
of £ can be identified with (—£, e) x £ in such a way that £ corresponds to {0} x £, that the 
metric takes the form dt 2 — gs where gs is independent of t and that V E is the pull-back of 
its restriction to £ under the projection (— e,s) x £ —> £. 



Fig. 1. Manifold X with product structure near £ 

In this case this representation is thought of as the preferred vacuum representation of an 
observer on £. 

3.6. Integral kernels of spectral projectors. As before let S© be the unique operator 
C^{X- S* R X <g> E*) -A C°°(X; S* L X <g> E*) such that 

o;s(4'(r’)'I , (a)) = f {v, S^u) dV 
Jx 

and se its integral kernel. Now let u G T 2 (£; S* R X ® E*). Then Gr^{u6y) G SoI(Vh) so that 
T(a) = B(Gr^{u5y) © 0) is defined. Similarly, \k(u) is defined for v G T 2 (£; SlX © E ). We 
compute 

ui£(*I'(i;)'I'(ii)) = ws(B(0 © Gr(v5^))B(Gr^(u5y,) © 0)) 


= ($7, 7Te(I?(0 © Gr(v5-£))B(Gr^(u5’e) © 0))f7)^^ 


= (ttz(B(Gr(v5z) © 0))17, tty,(B(Gr^(u5z) © 


(^Ppi:{Gr(v5z) © 0) A 17, Ppy,{Gr^{u5y) © 0) A 17 


HH S ) 


(8) 




ni 


= f (v,P>(%)^u))dA. 

J E 
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In other words, the integral kernel of the projector P>(Ve) coincides with s E (y, x)t^(x), 
restricted to (y, i) £ E x E, 

More generally, let Si and S 2 be two spacelike smooth Cauchy hypersurfaces and let Csi,e 2 
the unitary evolution operator defined in (5). For u £ L 2 ( Si; S^X <g> E*) we have by (8) that 
Gr,*((i^ 2 U^^ 2 u)Se 2 ) = and hence ^ r (^E 2 C/ E 1 1 ,E 2 #Ei' iX ) = 'M'M and similarly 

for \k. If v £ S 2 (Si; SlX <gj E ) we have on the one hand 

we 2 W^M) = / {v{y),ST, 2 {y,x)u{x))dX{x)dk(y) 

JEi xEi 

and on the other hand 


= (^-E 2 ^E 2 ^Ei 1 ,E 2 ? /Ei ?; ’-P>(^S 2 )( ? /E 2 ^E 2 ^ r Ei 1 ,E 2 ^Ei lt )) i2( . S2) 

= (^ r Ei 1 ,E 2 ^Ei' U ’P>(^S 2 )(^ r Si 1 ,E 2 ^Ei U ))L2(y; 2 ) 

= (4il J FE 1 ,E 2 P>(ys 2 )(^;E 2 ^ 1 «))L 2( E l) 

= [ {v,U^ 2 p>{^ 2 )U^Ji^u)dA. 

J Si 

Hence the integral kernel of Cei,e 2 °P> (^ 7 e 2 )°^Ei 1 e 2 coincides with ss 2 (y, x)if,^ ( x ), restricted 
to (y, x) £ Si x Si. 


3.7. Relative charge between Fock states associated to Cauchy hypersurfaces. Now 

we consider the following situation: let Si and S 2 be two spacelike smooth Cauchy hypersur¬ 
faces with Si lying in the past of S 2 . We assume that the metric of X and the connection 
V E of E have product structure near both hypersurfaces. We compute the relative charge 
for the Fock states ws 1 and we 2 : 


qW Si ,u>s 2 = / u El ( nSi ) dA 

J El 

= / tr (s El (a;,a;)-ss 2 (a;,a;))) dA(x) 

J Ei 

= / tr ((s Sl (x,a;) - s S2 (x,x))^^ ) dA(x) 

J Si 

= TV(p>(^ El ) - U^ 2 p>(n 2 )U^ 2 ) . 

It was shown in [7, Thm. 6.5] that the operator P>(^Ei)“^ r s 1 ,E 2 P>(^s 2 )^ r E 1 1 E 2 has a smooth 
integral kernel. In particular, the operator is of trace class. It now follows from [4, Thm. 4.1] 
that the trace is an integer and equals an index, namely putting U = t/ E] >E2 and pj = p>($7 E ), 

= ind [UP2U~ 1 P1 : PiS 2 (Si) -A Up 2 U~ 1 L 2 (E 1 )\ 

= ind [p 2 U- l Pl : piL 2 (Si) -a p 2 U~ 1 L 2 (E 1 )\ . 

By Theorem 4.1 and the concluding remark in [7] this index is given by 


(16) Q^ 1,WE2 = - / A A ch^) + 


My si)+My s 2 )+Mysi) - Mys 2 ) 


' M 


2 


h(y S2 ) 
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(17) 


= - / A A ch(V £ ) + 


'M 


~ ~ 7 ?(V , s 2 ) 

2 


Here M is the region between Si and £ 2 , he. M = J + (£i) D J (£ 2 ) where J + and J 
denote the causal future and past, respectively. By A we denote the A-form computed from 
the curvature of X and ch(V E ) is the Chern-character form for the curvature of V E , see [9, 
Sec. 4.1]. Hence A contains the contribution of gravitation to the relative charge and ch(V E ) 
that of the external field. 

Moreover, 77 (^ 2 ) denotes the 77 -invariant of the Dirac operator on the Cauchy hypersurface 
£ and hffix) the dimension of its kernel. The additional term -h(^-£ 2 ) in (16) is caused by 
a different convention for the spectral projectors in [7] concerning the eigenvalue 0. There is 
no trangression boundary term because of the product-structure assumption near £1 and £ 2 . 


3.8. Summary and examples. We summarize the results we have obtained. 


Theorem 3 . 1 . Let X be an even-dimensional globally hyperbolic Lorentzian spin manifold, 
let £i,£2 C X be two spacelike smooth Cauchy hypersurfaces with £1 lying in the past o/£ 2 . 
We put M := J+(£i) C J"(£ 2 ). 

Let E —> X be a Hermitian vector bundle with compatible connection X E . We assume that 
the metric of X and X E have product structure near £1 and £ 2 . 

Then the relative right-handed charge Q R x ’ as defined in (14) for the Fock states ujej 
and the Dirac operator twisted with E is given by 


(18) 


q ^ 2 



A A ch(V £ ) + 


~ h(y Sz ) + T](ff7 sj - r/ffigz) 

2 


Interchanging left-handed and right-handed spinors in the whole discussion we can also 
define the relative left-handed charge Q ^ 1 ' 2 • The projector defining the Fock representation 

is now given by P = p>(Vs) ©p>(Vs). Then a discussion analogous to the above yields 


(19) 


Q ^ 1 ’" 22 



A A ch(V E ) + 


-hCft s x ) + M^s 2 ) ~ 1) + ??(^s 2 ) 

2 


The exchange of chirality is equivalent to reversing the orientation of X. This explains the 
opposite sign for the contribution given by the integral. The induced orientations on £1 
and £ 2 will also be reversed which results in a replacement of ^y ij by — . Hence the 

77 -invariants get the opposite sign. The different convention concerning the eigenvalue 0 in 
the definition of the projection P is responsible for the opposite signs in the /t-terms. 

If we call i ^ 2 ;= 1^-2 the re i a n ve total charge and ■= 

Q U r1 ’ UJt “ 2 — Q L ^ Sl ’ UJS2 the relative chiral charge then (18) and (19) imply 


QWE^WEa = 0, 

QS ,W=2 = -2 [ A A ch(V' E ) + - h(y Sa ) + 77 (y Sl ) - 7 ?(^ 2 ). 

JM 

That is to say that the total charge is preserved in quantum field theory while the chiral 
charge is not. The chiral anomaly depends on the space-time curvature via A, on the external 
field via ch(V E ), and on spectral properties of the spatial Dirac operators via the h- and 
77 -terms. 
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Example 3.2. The following example is well known but shows that the contributions of 
the 77 -invariant are essential to give a correct integer valued total charge. Let I = M x S 1 
equipped with the metric dt 2 — dd 2 where 6 denotes the standard coordinate on the circle 
S 1 = R/LZ of length L. We twist with the topologically trivial complex line bundle E (indeed 
all complex line bundles on X are topologically trivial) and the connection V E = d+iA where 
the electro-magnetic potential A is of the form A = Ai{t)d9. The (real) curvature form of V E 
is given by F = dA = A\dt/\d9. This describes an electric held with no magnetic component. 

The surface X with the given metric has vanishing curvature. But even if it had nontriv¬ 
ial curvature it would not enter the formula for the relative charge because A has nonzero 
contributions only in dimensions divisible by 4. Thus A = 1. In two dimensions the Chern 
character form is simply given by ch(V E ) = F = ^dt A dO. 

For f 1 < t 2 and Ej = {tj} x S 1 we have M = [ii,^ 2 ] x S ' 1 and the integral is given by 

I A A ch(V B ) = f f p-d6dt= ^-(Ai(t 2 ) - Ai(ti)) . 

Jm Jt x Js 1 Z7T Z7T 

Now S 1 (and hence X) has two inequivalent spin structures, the trivial and the nontrivial 
spin structure. Spinors with respect to the trivial spin structure correspond to complex-valued 
functions. Spinors on S 1 with respect to the nontrivial spin structure correspond to functions 
u on R which are antiperiodic with period L , i.e. u(9 + L) = —u{6 ). The twisted Dirac 
operator on E = {t} x S l takes the form Vs = ^0 ~ Ai(t). Its eigenvalues have multiplicity 1 
and are given by 

27T . / \ 

—k - Ai(t) 

for the trivial spin structure and 


2 vr 

T 



Ai(t) 


for the nontrivial spin structure where k E Z. A little computation using Hurwitz ^-functions 
shows for the trivial spin structure 


2 


L 

2tt 


Ai(t) + h(rfx) + ??(y s ) = 2 


L 

2n 


Ai(t) 


+ 1 


and hence 






- 2 




For the nontrivial spin structure we obtain 


U2 Sl ,W£ 2 

^chir 


= 2 


L 

2tt 


A\(ti) 


1 

2 


-2 


L 

2tt 


Al(t2) 


1 

2 


Example 3.3. Let X = R x 5" 4fc_1 with a metric of the form dt 2 — gt where gt is a one- 
parameter family of Riemannian metrics on S’ 4fc_1 . The manifold X has a unique spin struc¬ 
ture. This time we let E be the trivial line bundle with trivial connection X E so that there 
is no external field. 

We consider M = [ti, t^\ x5 4fc_1 . It is shown in [7, Sec. 5] that the family g t and the Cauchy 
hypersurfaces Ej = {tj} x 5' 4fc_1 can be chosen in such a way that 1 ) = M^e 2 ) = ^ an< i 
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f M a + ^A1 = (_!)« ffi. Hence 

<?£;■"” = (- d‘2(“). 

In this case the space-time metric, i.e. gravity, is causing the nontriviality of the chiral anom¬ 
aly. In dimension 4 {k = 1) this example was already discussed in [22], 

Example 3.4. Next we consider Bianchi-Type-I space-times. Here I = lx R 3 /r where T is 
a lattice in R 3 . Denoting the coordinate of the R- factor by t and the standard coordinates of 
R 3 by x 3 , a? 2 , x 3 the metric of X takes the form dt 2 — a 2 (t)(dx 1 ) 2 — a\(t)(dx 2 ) 2 — a 2 (t)(dx 3 ) 2 
where : R —> R are smooth positive functions. Hence, for fixed t, the Cauchy hypersurface 
{i} x M 3 /r is a flat torus. 

The A-integrand can be computed to vanish identically on X, A = 0. The manifold X 
has 8 different spin structures. The Dirac spectrum of a flat torus is symmetric about 0 for 
every choice of spin spin structure so that all ^-invariants vanish. Moreover, the kernel of the 
Dirac operator on a flat torus is 1-dimensional for one of the spin structures and trivial for 
the others. Thus the ^contributions in the formulas for the chiral anomaly cancel in 

all cases. We conclude that for Bianchi-Type-I space-times the chiral anomaly vanishes, 

Q^ 1,WEa = Q^ 1 ’^ 2 = QS ,WEa = °- 


Example 3.5. Finally we consider Bianchi-Type-II space-times. The torus in a Bianchi-Type- 
I space-time is replaced by a 3-dimensional Heisenberg manifold, i.e., X = R x (He(Z)\He(M)) 

/1 x z\ 

where He(R) and He(Z) denotes the Heisenberg group of matrices of the form I 0 1 y I 

V° 0 V 

with real or integral entries, respectively. Clearly, He(R) is diffeomorphic to R 3 but the 
quotient X := He(Z)\He(R) is a compact 3-manifold non-diffeomorphic to a torus. It has 
four different spin structures. For positive numbers a and b we obtain a left-invariant metric 
on He(R) by 


9a,b 




dx 2 + 



+ a 2 ^j dy 2 + b 2 dz 2 


b 2 xy 


dxdy + b 2 ydxdz — b 2 xdydz . 


By left-invariance the metric descends to a Riemannian metric on the quotient X. 

Now we let a, b : R —>• R be smooth positive functions and we equip X with the metric 
dt 2 — g a (t),b(t )• ^ tedious but straightforward computation yields 


- {a 2 ba 2 — a 3 ba — a 3 ab + a 4 b — b 3 )(ba 

- = 48^5 

Integration over M = [t\,t 2 ] x X yields 


ab) 


dt A dx A dy A dz . 



2 a 2 b 2 a 2 + 4 a 3 abb 
192vr 2 a 4 


2 a A b 2 


t2 


t=t\ 


Assuming product structure near t± and t^ as required in Theorem 3.1, the derivatives of a 
and b vanish at t\ and at t 2 so that we are left with 


A = 


1 f b 4 (t i) 

1927T 2 \a 4 (ti) 


b 4 (t 2 ) \ 

a 4 (t 2 )J 


M 
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The ? 7 -invariant of the Dirac operator on is computed in [25, Thm. 11] (with A = a 2 /b 2 , 
r = vi = W 2 = rn v = 1 and u 2 = w\ = m w = 0 in their notation): 

^ = WT)- m 

where N(t ) is an explicitly given integer depending on a(t)/b(t) and the spin structure. From 
the explicit formula (44) in [25] one sees that for every choice of spin structure N(t) —>• oo if 
b{t)/a{t) —>• oo and N(t) —> 0 if b(t)/a(t ) —> 0. Therefore the functions a(t) and b(t ) can be 
chosen such that IV(i 2 ) — N(t\) is arbitrarily large. Moreover, for generic choice of a(tj)/b(tj) 
we have = 0. Then 

Qchir’^ 2 =-2 [ & + h(y Sl )-h(f^) + r,(t Sl )-r,(f Sa ) = N(t 2 )-N(t 1 ) 

Jm 

can take arbitrarily large positive but also arbitrarily small negative values for suitable choices 
of the functions a(t) and b(t). 


4. Conclusion and outlook 

Given two Hadamard states on the field algebra for the Dirac equation, we defined a relative 
current which can be thought of as the expectation value (with respect to one state) of the 
normally ordered current operator where normal ordering has been done with respect to the 
other state. In contrast to currents associated with one state, these relative currents can be 
defined unambigously without any need for regularization. Integration of the relative current 
over a Cauchy hypersurface yields the relative charge. In case the two states are vacuum 
states of observers moving along two distinguished timelike Killing fields, we expressed the 
relative charge in terms of an integral into which the gravitational field and the gauge field 
enter plus a boundary contribution essentially given by the r/-invariant. The crucial tool was 
a novel index theorem for hyperbolic Dirac operators which makes it possible to avoid any 
Wick rotation thus making the whole derivation mathematically rigorous. 

Of course, the assumption on the space-time to be spatially compact is something that one 
would like to get rid of in order to better understand the local properties of anomalies and 
in order to be able to treat more models from general relativity. Technically, this assumption 
ensures that the relative charge is finite and can be expressed as an index. However, the 
relative current is locally defined and does make sense also on more general space-times. A 
local index theorem should be able to compute the local current generated by the external 
field. In case of noncompact Cauchy hypersurfaces which have compact quotients there is 
also the possibility to employ the concepts of L 2 -indices, L 2 -traces, and L 2 -dimensions as 
introduced by von Neumann. 

Appendix A. The canonical trace on the algebra of smoothing operators 

WITH COMPACT SUPPORT IN ONE VARIABLE 

Let X be a manifold with a positive volume density dV. Suppose that E is a complex 
vector bundle over X. Then any smooth section a of the bundle EM E* will determine an 
operator A : Cfi°(X;E) -A- C°°(X;E ) via 
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If the integral kernel a has support contained in a set of the form K x X, where K C X is 
compact we say that a is compactly supported in the first variable. The set of operators with 
integral kernel that is compactly supported in the first variable shall be denoted by T^ 00 . 
Note that this set forms an algebra and these operators map Cq°(X: E) to Cq°(X', E). A 
canonical trace can be defined on this algebra by 

Tr(A) := / tr a(x, x) dV(x). 

J x 

This trace cannot be directly interpreted as an L 2 -trace because an operator in Tq 00 need not 
necessarily be L 2 -bounded. However, given a compact exhaustion K n of X and a sequence 
of compactly supported smooth functions Xn with Xn( x ) = 1 for all x G K n , we have 

Tr(A) = lim Tr i2 (Ax n ). 

n—> 00 

The operator Ax n has smooth integral kernel a(x,y)x n (y ) which is compactly supported in 
both variables. Hence it is trace class by Mercer’s theorem. For A, B G Tg 00 the trace 
property Tr (AB) = Tr (BA) follows immediately from Fubini’s theorem. 

Now let P be a differential operator acting on sections of E. If A G d'g" 00 then PA, AP G 
Tg°° because they have integral kernels (P <8>1e*)ci and (l#® P*)a, respectively. The trace 
property Tr(PA) = Tr(AP) also holds even though P ^ 'f'o" 00 - This can easily be checked 
using integration by parts and the fact that in local coordinates 

d d d 

d^ a ^ x,x) = d? a{ - x ' y) \ y=x + fr/ a ( x ' y ^y= x - 
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